Transversely Stable Soliton Trains in Photonic Lattices 
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We report the existence of transversely stable soliton trains in optics. These stable soliton trains 
are found in two-dimensional square photonic lattices when they bifurcate from X-symmetry points 
with saddle-shaped diffraction inside the first Bloch band and their amplitudes are above a certain 
threshold. We also show that soliton trains with low amplitudes or bifurcated from edges of the 
first Bloch band (F and M points) still suffer transverse instability. These results are obtained in 
the continuous lattice model and further corroborated by the discrete model. 
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Oh' It is well known that in homogeneous optical media, a bright soliton stripe, which is localized along the longitudinal 
' direction and uniform along the transverse (stripe) direction, is always unstable to transverse perturbations [iHtI- 
When the longitudinal and transverse diffractions have the same sign, the instability is of neck-type which breaks up 
T— I ' the stripe into filaments [1, S 01 j while when the two diffractions have the opposite sign, the instability is of snake-type 
, which bends the stripe to a wavy form 0-0|- These two types of instabilities have been experimentally observed as 

■ well [3, 8]. When a one-dimensional optical lattice is imposed along the longitudinal or transverse direction of the 
' soliton stripe, the soliton stripe is still transversely unstable [1, While the transverse instability can be utilized 
' certain applications (such as pulse compression), in many other cases (such as experiments in lower dimensions), it 

is detrimental and undesirable. To suppress this transverse instability, some ideas have been proposed. For instance, 
this instability can be completely eliminated if the soliton stripe is made sufficiently incoherent along the transverse 
direction Tlj. This instability can also be significantly reduced (but not eliminated) by nonlinearity saturation 
^ or incoherent mode coupling jl2l. [isj. Beside optics, transverse instability is also a common phenomenon in other 
branches of physics such as water waves [1, [l^-[l6| . 

In this paper, we report the existence of coherent and transversely stable solitons trains in optics (for the first time 
^ to our knowledge). These stable soliton trains are discovered in two-dimensional photonic lattices. They comprise a 

O ^ periodic sequence of intensity lumps along the transverse direction, and are localized along the longitudinal direction. 
These soliton trains, which are exact stationary solutions of the underlying nonlinear lattice system [l7l - [l9| . are 
transversely stable when they bifurcate from the X-symmetry point (with saddle-shaped diffraction) inside the first 
Bloch band of the lattice, and the amplitude of the soliton trains is above a certain threshold value. This finding is 
Ijr^ ' first obtained in the full continuous lattice model. Then it is also corroborated by the discrete nonlinear Schrodinger 
model, which is shown to support transversely stable discrete line solitons as well under similar conditions. Physically, 
\^ •- these soliton trains are transversely stable due to the stabilizing effect of the photonic-lattice potential. However, the 
I photonic lattice does not stabilize every soliton train. Specifically, low-amplitude soliton trains as well as soliton trains 

■ that bifurcate from edges of the first Bloch band (F and M points) are still transversely unstable. These transversely 
unstable soliton trains will also be explained both mathematically and physically. 

II. SOLITON TRAINS IN TWO-DIMENSIONAL PHOTONIC LATTICES 

, The theoretical model we use for coherent-beam propagation in a two-dimensional (2D) photonic lattice is Q 

iU, + U^^ + Uyy + n{x,y)U + a\U\'U = 0, (2.1) 

where z is the direction of propagation, {x, y) is the plane orthogonal to the propagation direction, n{x, y) is the 
periodic refractive-index variation on the orthogonal plane, and a = ±1 represents self-focusing and self-defocusing 
nonlinearity. All variables have been normalized. In our analysis, we take the lattice to be 

n(x, y) — h{cos^x + cos^y), (2.2) 

where h is the index-variation depth parameter. This is a square lattice which arises frequently in optics and Bose- 
Einstein condensates 7,, 20]. This lattice is 7r-periodic along both x and y directions, and is displayed in Fig. [ija). 
Below, X and y will be called the principal axes of the lattice since along them the lattice has the smallest period. 
This lattice supports soliton trains aligned along various directions in the {x,y) plane pTt - fToj . In this paper, we only 
consider soliton trains alighed along a principal axis of the lattice for simplicity, and take this principal axis to be the 
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FIG. 1: (Color online) (a) the photonic lattice (|2.2p (with h = 6); (b) the diffraction surface fx = fi{ki, ^2) of the first Bloch 
band for this lattice. The high-symmetry points F = (0,0), X = (1,0) and M = (1,1) in the irreducible Brillouin zone 
— 1 < fci, fc2 < 1 are marked by red dots. 



y-direction. Such soliton trains are of the form 

Uix,y,z) = uix,y)e-'''\ (2.3) 

where u{x, y) is a real-valued function which is localized along the longitudinal x-direction and periodic along the 
transverse y-direction, and /i is the propagation constant. The function u{x,y) satisfies the equation 

Uxx + Uyy + fiu + n{x, y)u + au^ = 0. (2-4) 

In this section, we examine these soliton train solutions. 

When u{x,y) is infinitesimal, Eq. (|2.4p becomes a linear equation whose bounded solutions are Bloch modes, 

p{x,y;fi) = e'^''^-+''^y'>p{x,y;fi), (2.5) 

where p(x, y; /i) is 7r-periodic in both (x, y), ki, k2 are wavenumbers in the irreducible Brillouin zone —1 < fci, k2 < 1, 
and 

/i = ^(fci,fc2) (2.6) 

is the diffraction relation. At the lattice depth h = 6, this diffraction function for the first (lowest) Bloch band is 
shown in Fig. [Hb). If the linear Bloch mode p(x, y; /i) is required to be real, this mode then must be at one of the 
high-symmetry points, F = (0,0), M = (1, 1), and X — (1,0), (0, 1) in the irreducible Brillouin zone, and be tt- or 
27r-periodic in x and y. For the first Bloch band (see Fig. mb)), F and M are the lower and upper edge points, and 
X lies inside the Bloch band where the diffraction surface has a saddle shape. In this paper, we only consider soliton 
trains bifurcated from this first band. 

When u{x,y) is small but not infinitesimal, Eq. (j2.4p is weakly nonlinear, and its solution is a weakly modulated 
Bloch-wave packet which bifurcates out from the underlying high-symmetry point of the Bloch band. When this 
solution is a soliton train along the y direction, it can be expanded into a perturbation series 

uix, y) = eA{X)p{x, y) + e''A'{Xy^ (x, y) + . . . , (2.7) 



/i = ^in + re^ (2.8) 

where ^0 is the propagation constant of the high-symmetry point, p(x,y) is the Bloch wave at fiQ, vi{x,y) is a 
generalized Bloch function at /iq which satisfies the equation 

dxxVi + dyyVi + [/io -I- n{x, y)]vi = -2dxP, (2.9) 

T = ±1, < e ^ 1, X = ex is the slow spatial variable, and A{X) is the one-dimensional (ID) envelope of this Bloch 
wave. Notice that the solution vi{x,y) to Eq. (|2.9p is periodic in both x and y with the same period as p{x,y). 
In addition, this solution is not unique since one may add an arbitrary homogeneous solution ^p(x,y), where C is a 
free constant. Returning to the expansion (j2.7p . we can see that adding C,p{x,y) to vi^x^y) amounts to a shift in 
the position of the envelope A{X). In order to fix the location of the envelope A{X), we require that vi{x,y) be 
orthogonal to p{x, y), 

p{x,y)i'i{x,y)dxdy ^0. (2-10) 
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This orthogonality requirement uniquely determines the solution Vi{x^y). In the present case, the lattice n{x,y) in 
(|2.2p is symmetric in x, then p{x,y) is either symmetric or antisymmetric in x. In this case, ui{x,y) would have the 
opposite x-symmetry of p(a;, y) under the above orthogonality condition. 

By inserting this perturbation expansion into Eq. (j2.4p and following the analysis very similar to [2l| , we find that 
the envelope A{X) satisfies the following equation 

DiAxx+rA + aaA^ =Q, (2.11) 

where the x-direction diffraction coefficient Di and the nonlinear coefficient a are given by 

j\ ^Z^p'^{x,y)dxdy 

ct = . 

M=Mo' Il^Il^P^ix,y)dxdy' 



_ l 9V(fcl,fc2) 
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When sgn(_Di) = sgn{a) — — sgn(T), the envelope equation (j2.11|) admits a sech soliton 



A{X) = . / ^ sech ^^° , (2.12) 

where Xq — exo is the location of the peak of the envelope function A{X). When the Bloch wave p{x, y) is modulated 
by this ID sech envelope, the resulting solution (|2.7|) is then a low-amplitude soliton train along the y direction. 

One may notice that the envelope equation (|2.1ip is translation-invariant, hence Xq is a free parameter in the 
envelope solution (j2.12p . This seems to imply that soliton trains can be obtained regardless of the position of the 
envelope (12.121) relative to the underlying periodic potential. This is not true however. In a similar ID-lattice model, 
it has been shown that the peak of the envelope can only be located at two positions relative to the lattice 0, [H, [l^ . 
Slight extension of that analysis to the present 2D soliton train problem shows that the envelope p.l2p of the soliton 
train, for the symmetric cosine lattice (|2.2p . also must be located at one of the two positions 



Xq = 0, 7r/2. (2.13) 

The resulting soliton train with a;o = is called the on-site soliton train, and the other one with xq = 7r/2 is called the 
off-site soliton train. The off-site train resides at (x, y) regions of low refractive indices n{x, y) and is thus expected 
to be always unstable 0, [H, H^. Thus we only consider on-site soliton trains in the rest of this article. 

When the amplitude of the soliton train is not small, the above perturbation series would be invalid. In such cases, 
soliton trains can be determined numerically ff\. To illustrate, we take self-defocusing nonlinearity (ct — —1) and 
lattice depth h — Q. In this case, a family of soliton trains bifurcates out from the X point inside the first Bloch 
band. The power curve of this family is displayed in Fig. HJa) . Here the power P is defined as the integral of v? from 
—oo < X < oo and on one transverse period < ?/ < tt, i.e., 

/OO /"TT 
dx I dy u^{x,y). (2.14) 
-oo JO 

Two typical soliton trains, with low and high amplitudes (powers), are displayed in Fig. [2jb,c) respectively. The 
propagation constants for these two soliton trains are fi = 4.3495 and 5.6. One can see that the intensity peaks of 
these solitons are in-phase along the transverse {y) direction and out-phase along the longitudinal (x) direction. In 
addition, the soliton train near the X point has low amplitude and power, and is longitudinally broad (occupying 
many lattice sites), while that away from the X point has high amplitude and power, and is longitudinally strongly 
localized (occupying practically a single lattice site). This type of soliton trains has been theoretically predicted 
and experimentally observed in jl9l] . In Fig. |4] of later text, another family of soliton trains bifurcated from the 
M-symmetry point of the first Bloch band will also be displayed. 

In the next section, we will examine the transverse stability of soliton trains. Before detailed analysis, let us first 
develop some intuition. The soliton trains exist under both self- focusing and self-defocusing nonlinearities (cr = ±1), 
and they can bifurcate from the F and M points at edges of the Bloch band [13, [3 , or from X points inside the 
Bloch band [l9'| . For soliton trains bifurcated from the F point (which occurs under self- focusing nonlinearity) [it'] , the 
intensity peaks along the transverse (train) direction are all in-phase. For soliton trains bifurcated from the M point 
(which occurs under self-defocusing nonlinearity) p^], the adjacent intensity peaks along the transverse direction are 
all out of phase (if this transverse direction is along a principal axis of the lattice as in our present case, see Fig. 
3]). It is known that in-phase dipoles under self- focusing nonlinearit y a nd out-phase dipoles along a principal axis 
of the lattice under self-defocusing nonlinearity are both unstable [tI. |24| - [27| . Then if the soliton trains are strongly 
localized along the longitudinal direction (which occurs at large amplitudes, see Fig. |H[c)), we may view the soliton 
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FIG. 2; (Color online) (a) the power curve of soliton trains bifurcated from the X point inside the first Bloch band under 
defocusing nonlinearity; the dashed (red) segment is transversely unstable, while the solid (blue) segment is stable; (b,c) profiles 
u{x, y) of soliton trains at low and high amplitudes respectively; these solutions are located at the points marked by the same 
letters on the power curve of (a); (B,C) linear-stability spectra for the soliton trains in (b,c). 

trains as a collection of transverse dipoles, hence we may expect the soliton trains bifurcated from V and M points 
to be unstable. However, for soliton trains bifurcated from the X points (under either self-focusing or self-defocusing 
nonlinearity), adjacent transverse intensity peaks are out of phase under self- focusing nonlinearity, and are in-phase 
under self-defocusing nonlinearity (see Fig. [5] and fl9'|). Dipoles with such phase structures are stable in deep lattices 
^, 24-27]. Thus soliton trains bilFurcated from X points (at high amplitudes as in Fig. ^c)) may be free of transverse 
instabilities. In the next two sections, we will confirm that these intuitions are largely correct, hence transversely 
stable soliton trains will be discovered. 

It should be cautioned, however, that these intuitions are reasonable only when the soliton trains have high ampli- 
tudes, so that they are strongly longitudinally localized (as in Figs, and|H[c)), which makes the transverse-dipole 
analogy meaningful. If the soliton trains have low amplitude, then they are longitudinally broad and occupying 
many lattice sites (see Fig. [DJb)). In that case, longitudinal inter-site coupling becomes important, which makes the 
above transverse-dipole analogy inappropriate. As we will show later, all low-amplitude soliton trains are transversely 
unstable. 



III. TRANSVERSELY STABLE SOLITON TRAINS 



In this section, we study the transverse stability of soliton trains in 2D photonic lattices, and show that transversely 
stable soliton trains exist. 

First, we briefly show that low- amplitude soliton trains are always transversely unstable. To show this, we consider 
the dynamics of a low-amplitude soliton train, whose solution can be expanded into a perturbation series, 

U{x, y, z) = e-*'^"^ [e*(X, Y, Z)p{x, y) + e^U2 + . . .] , (3.1) 

where /iq is the propagation constant of a high-symmetry point, p{x,y) is the Bloch wave at /io, < e <ti 1, X = ex, 
Y = ey are slow spatial variables, Z = e^z is the slow propagation-distance variable, and 'i>{X,Y, Z) is the envelope 
function of this low-amplitude soliton train. Following the analysis of [2l|, it is easy to show that the evolution of the 
envelope function ^(X, Y, Z) is governed by the following constant-coefhcient 2D NLS equation 

z^fz + i'l^'xjf + £'2l'Yy + cral^'p^' = 0, (3.2) 

where D2 is the y-direction diffraction coefhcient. 
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and Di , a have been given before. This 2D envelope equation admits a hne-sohton solution 

'f{X,Y,Z) = A{X)e-"^, (3.3) 

where r = —a, and A{X) is given in Eq. (|2.12|) . This envelope line soliton, when modulated onto the Bloch mode 
p{x, y), yields the low-amplitude soliton train (j2.7l) derived in the previous section. It is well known that this envelope 
line soliton (|3.3I) is transversely unstable in the constant-coefScient 2D envelope equation (13. 2|) [H, 0, • Thus low- 
amplitude soliton trains (|2.7|) are also transversely unstable. Since the transverse instability of low-amplitude soliton 
trains is induced by the transverse instability of their envelope line solitons, obviously the longitudinal wave coupling 
plays an important role in this instability (as has been pointed out in the end of the previous section) . 

We can further derive the analytical formula for unstable eigenvalues of low-amplitude soliton trains. Since the 
soliton train is periodic along the transverse direction, the normal modes of infinitesimal disturbances to the train are 
of the form ^ ^ e^^y'^^^ip{x,y), where k is the transverse wavenumber of the disturbance, A is the eigenvalue, and 
ip^x, y) is vr-periodic in y. From the above analysis and after variable scalings, it is easy to show that the eigenvalue 
A(fc) is given by 

X{k)=e'Ai^\D^\k/e), (3.4) 
where A{K) is the eigenvalue of the normalized sech line soliton 

$o(^, Y,Z)^V2 sechXe'^ (3.5) 

in the unit- coefficient 2D NLS equation 

i<i>z + + sgn{DiD2)^YY + l^'P* = (3.6) 

for the disturbance proportional to e*^^, which have been obtained in [3, [1^. A more rigorous derivation of the 
eigenvalue formula p.4[) . which is based on the study of the linear-stability eigenvalue problem of low-amplitude 
soliton trains, can also be made by a slight modification of the analysis in [29|. 

If the low-amplitude soliton train bifurcates from the X-symmetry point (as in Fig. [2|), sgn{DiD2) = -1, hence the 
envelope's transverse instability is of snake type. Since D1D2 < 0, the diffraction surface at the X point has a saddle 
shape (see Fig. [Ifb)). If the train bifurcates from the edges of the first Bloch band (as in Fig. sgn(£>il?2) = 1, 
hence the envelope's transverse instability is of neck type. 

When the amplitude of the soliton train is not small, the analytical stability analysis above becomes invalid. 
Below, we numerically determine the transverse stability of higher-amplitude soliton trains (using numerical methods 
described in [t'I). 

First, we consider the soliton trains in Fig. [2l which bifurcate from the X-symmetry point of the first Bloch band. 
The numerically obtained linear-stability spectra for the low- and high-amplitude soliton trains of Fig. [2jb,c) are 
displayed in Fig. [2][B,C) respectively. The spectrum in Fig. [2]^B) contains both real and complex unstable eigenvalues 
which lie on the right half of the spectral A-plane, indicating that the low-amplitude soliton train in Fig. [H^b) is 
transversely unstable. This numerical result agrees with the analytical result given above. Quantitatively we have 
also compared these numerical eigenvalues with the analytical formula p.4|) . with the function A{K) provided in 
0, [13 , and found excellent quantitative agreement as well. 

A more important finding on this family of soliton trains is that, when the amplitude (or power) of these soliton 
trains reaches above a certain threshold, transverse instability disappears, and these soliton trains become fully stable. 
This can be seen in Fig. [2{C), which gives the linear-stability spectrum for the high-amplitude soliton train in Fig. 
m^c). This spectrum does not contain any unstable eigenvalue, indicating that this high-amplitude soliton train is 
fully stable. What happens is that, as fj, moves away from the X point /xq — 4.330 (i.e., as e increases), the unstable 
eigenvalues initially increase in size and drift away from the imaginary axis of the spectral plane, as is predicted by 
the analytical formula p.4|) . However, as fi increases further, the unstable eigenvalues turn around and start to move 
toward the imaginary axis (hence the transverse instability weakens). When ^ > ~ 5.44, or P > 5.3, all unstable 
eigenvalues merge into the imaginary axis, hence instability vanishes, and the soliton trains become transversely 
stable. Physically, what happens is that, as ^ moves away from the X point, the soliton train transforms from a 
low-amplitude longitudinally-broad train into a high-amplitude longitudinally- narrow train (see Fig. [ll^b,c)). In this 
process, the snaking instability of the low-amplitude longitudinally-broad soliton train is eventually arrested by the 
lattice when the train becomes strongly localized longitudinally. 

We have found that the high-amplitude soliton train in Fig. (^fc) is not only linearly stable, but also nonlinearly 
stable. To demonstrate its nonlinear stability, we simulate its evolution in Eq. (j2.ip when it is perturbed by 10% 
random-noise transverse perturbations. After long-distance simulations, we have found that this train remains robust 
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FIG. 3: (Color online) Nonlinear evolution of the soliton train in Fig. [2]^c) under 10% initial transverse perturbations (shown 
are intensity fields |f7p). Left: the initial perturbed soliton train; right: the solution at z — 100. 



and does not break up at all. For instance, one simulation result (with simulation distance z = 100) is displayed in 
Fig.H 

The soliton train in Fig. ^c) may remind us the train of lumps which form after the onset of neck-type transverse 
instability to a line soliton in a homogeneous medium [3,0|- However, it is important to recognize that the train of 
lumps in a homogeneous medium is not a stable structure. Upon further propagation, it will break up Q- In contrast, 
the soliton train in Fig. [2jc) is both linearly and nonlinearly stable, and can propagate for all distances without 
breakup (see Fig. [3]). Here the photonic lattice plays an important role in the stabilization of the soliton train in the 
presence of transverse perturbations. 

The photonic lattice, however, cannot stabilize every soliton train (even if their amplitude is high). To demonstrate, 
we consider the soliton trains which bifurcate from the upper edge of the first Bloch band (M-symmetry point) under 
self-defocusing nonlinearity (with ft, = 6 as before). This soliton family is displayed in Fig. |ll^a,c). From the 
solution profile shown in Fig. Sfc), we see that the intensity peaks of these soliton trains are out of phase with 
each other along both the transverse and longitudinal directions (as is expected since these trains bifurcate from 
the M-symmetry point). However, this family of soliton trains are all transversely unstable. To demonstrate, the 
linear-stability spectrum for the high- amplitude soliton train in Fig. |3^c) is shown in Fig. Hljb). The real unstable 
eigenvalues on the right half of the spectral plane indicate that this high-amplitude soliton train is linearly unstable. 
The nonlinear instability of this soliton train is displayed in Fig. Hl^d). It is seen that under weak perturbations, 
this soliton train breaks up into filaments. For this solution family, the transverse instability of low-amplitude soliton 
trains is neck-type (since the diffraction coefficients Z?i,Z?2 have the same sign at the M point). Then Fig. 0] shows 
that the neck-type transverse instability of low-amplitude soliton trains cannot be arrested by the photonic lattice as 
the soliton train's amplitude becomes high. 

From Figs. [2]and|4l we learn that the phase relation of transverse intensity humps also plays an important role in 
the stabilization of soliton trains. Specifically, under defocusing/focusing nonlinearity, the transverse intensity peaks 
of the soliton train should be in-phase/out-phase in order for it to be stable (as one would expect from dipole stability 
results 0, [13-1131 )■ This requirement on the phase relation translates into a requirement that these soliton trains 
should bifurcate from the X-symmetry point inside the Bloch band. 

We should add that, for the family of soliton trains in Fig. [21 even though they become stable when /x > /Xc ~ 5.44, 
if ^ gets close to the second Bloch band (whose lower edge is at /x = 7.23), then these soliton trains would become 
unstable again due to wave coupling to the second band. 

IV. TRANSVERSELY STABLE LINE SOLITONS IN THE DISCRETE NLS EQUATION 

In the previous section, we demonstrated the existence of transversely stable soliton trains in the continuous lattice 
model (|2.ip . In this section, we corroborate this finding by showing that the analogous transversely stable line solitons 
exist in the discrete NLS model as well. The discrete NLS model is often used to qualitatively describe wave dynamics 
in the continuous lattice model (j2.ip with a deep lattice potential under the tight-binding approximation (27j . In this 
context, the variable in the discrete model can be viewed as the complex amplitude of the ground-state eigenmode 
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FIG. 4: (Color online) (a) the power curve of the soliton trains bifurcated from the M point of the first Bloch band under self- 
defocusing nonlinearity (the dashed red line indicates this whole family is transversely unstable); (b,c) the stability-spectrum 
and profile u{x, y) of the soliton train at the marked point on the power curve (/i = 5.6); (d) nonlinear evolution [z = 7) of the 
soliton train in (c) under 10% initial transverse perturbations (shown is the intensity field |f7|^). 



of each lattice-cell potential. The approximation of the continuous model (|2.ip by the discrete NLS equation is a 
significant reduction. Under this reduction, the soliton train shown in Fig. [Hb,c) in the continuous model becomes 
a transversely uniform line soliton in the discrete model (see Fig. [5jc,d) below). This treatment simplifies both 
analytical analysis and numerical computations. 
The discrete NLS equation we consider is 

i^Um,n + ^2Urn.n + Cr| i7„,„ ^ [/,„,„ = 0, (4.1) 

az 

where A2 is the two-dimensional difference operator 

(^m,n-t-l '^Ujyi,n ~i~ ^m.n— l); 

and CT = ±1 is the sign of nonlinearity. Here the intersite-coupling coefficient in front of lS.2Um,n has been normalized 
to be unity through scalings of z and Um,n- An important property of this discrete model is that self-focusing and 
self-defocusing nonlinearities can be transformed to each other since this model is invariant under the transformation 

M 

«7™,n ^ (-l)™+"e-8-C/;,,„, a^-cj. (4.2) 

This transformation is very helpful for us to understand the connection on soliton configurations and their stability 
properties between self-focusing and self-defocusing nonlinearities. It also means that one only needs to study one 
type of nonlinearity (say, self-defocusing nonlinearity), and infer the results for the other type of nonlinearity by this 
transformation. 

Solitons in the discrete model (|4.ip are sought in the form 

U^,n{z) = u™,„e-^''^ (4.3) 

where Um,n is a real-valued function which satisfies the equation 

^2Urn,n + M^m^n + CrMm,™ = 0, (4.4) 

and /i is the propagation constant. When Um,n is infinitesimal, the nonlinear term in ()4.4|) drops out, and the bounded 
solution to the remaining linear equation is then a discrete Fourier mode, Um,n = e*('^i'"+''2")^ where —it < ki, k2 < tt 
are wavenumbers along the m- and n-directions. Inserting this discrete Fourier mode into the linear part of Eq. (14. 4p . 
we get the diffraction relation 



fj. — 2(2 — cos fci — cos ^2). 



(4.5) 
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FIG. 5: (Color online) (a) the diffraction surface (|4.5p of the discrete NLS equation; (b) the power curve of discrete line 
solitons bifurcated from the X point inside the continuum band under self-defocusing nonlinearity; the dashed (red) segment 
is transversely unstable, while the solid (blue) segment is stable; (c,d) profiles Um,n of discrete line solitons at low and high 
amplitudes respectively; these solitons are located at the points marked by the same letters on the power curve of (b); (C,D) 
linear- stability spectra for the discrete line solitons in (c,d). 



The corresponding diffraction surface is shown in Fig. EJa). It is easy to see that this diffraction surface closely 
resembles the first Bloch band of the continuous model (see Fig. [1]). Thus the discrete model (|4.ip is appropriate 
for describing wave dynamics associated with the first Bloch band in the continuous model (|2.1[) . Notice that the 
continuous spectrum of this discrete model (|4.ip contains only a single band < /i < 8, while the continuous spectrum 
of the continuous model (12. ip often contains multiple Bloch bands. Thus if wave dynamics in the continuous model 
(12. ip involves higher Bloch bands, the discrete model (14.11) will be inappropriate. 

Corresponding to the continuous soliton trains in Fig. [21 we can find a family of discrete line solitons that bifurcate 
from the X-symmetry point of the continuum band under self-defocusing nonlinearity (cr ~ —1). These discrete line 
solitons Um.n are n-independcnt, thus they are simply ID discrete solitons in the 2D model (|4.4p . The power curve of 
this soliton family is shown in Fig. [5{b). Here the power P is defined as 

oo 

p{^i)^ E (4.6) 

m— — oo 

Profiles of two typical solitons, near and far away from the X point (with /i — 4.4 and 29), are displayed in Fig. 
[5Uc,d). The one in Fig. [Sfc) has low amplitude (power) and is longitudinally broad, which is the counterpart of the 
low- amplitude soliton train in Fig. [H^b). The one in Fig. [S{d) has high amplitude (power) and is longitudinally 
strongly localized, which is the counterpart of the high-amplitude soliton train in Fig. [2Jc). To determine the 
transverse stability of these discrete line solitons, we have computed their linear-stability spectra, and the results are 
displayed in Fig. [5jC,D). The spectrum in Fig. indicates that the low-amplitude discrete line soliton in Fig. 

[5]Jc) is transversely unstable. In addition, this spectrum qualitatively closely resembles the continuous counterpart 
in Fig. [IjB). The spectrum in Fig. [S^D), on the other hand, does not contain any unstable eigenvalue, indicating 
that the high-amplitude discrete line soliton in Fig. [Sfd) is transversely stable. Thus the existence of transversely 
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stable discrete line solitons is established. In the present discrete model (I4.ip . the threshold for transversely stable 
line solitons is at /i > /x^ ~ 21.5, or P > 39.0. 

Now we examine nonlinear developments of these discrete line solitons under transverse perturbations. We find 
that when these solitons are linearly unstable (see Fig. [S{b)), then under perturbations, they would develop snake 
instability and eventually disperse away. This is illustrated in Fig. [6l^a,b) for the low-amplitude discrete line soliton 
in Fig. [SJc). Analytically, this snake instability can also be understood. Briefly speaking, the envelope of a low- 
amplitude discrete solution to Eq. (|4.1I) is governed by an equation similar to p.2p . where Di and D2 are the 
diffraction coefficients. At the X-symmetry point, D1D2 < 0, thus the line soliton to this envelope equation suffers 
snake-type instability [ll,[l,|3|, which translates to the snake instability observed in Fig. |6l[a,b). When the discrete line 
solitons are linearly stable (see Fig. [SJb)), however, they would propagate robustly for all distances without breakup. 
This is illustrated in Fig. Hl^Cjd) for the high-amplitude discrete line soliton in Fig. [Sfd). 




-4 4 -4 4 

m m 



8 


z=0 




(c) 




z=100 




(d) 


n 

-8 


















-4 





4 


-4 





4 



m m 

FIG. 6: (Color online) Nonlinear evolutions of discrete line solitons in Fig. [5jc,d) under 10% initial transverse perturbations 
(shown are intensity fields |(7|^). Upper row: the low-amplitude case; lower row: the high-amplitude case. 

By comparing these results for the discrete NLS model (|4.ip to those for the continuous model (|2.1|) , one can easily 
see that the results for both models are qualitatively almost identical. For both models, we discovered transversely 
stable soliton trains under similar conditions, i.e., when they bifurcate from X-symmetry points and have high 
amplitudes. The instability behaviors for low-amplitude soliton trains are also the same in both models. However, 
minor differences between the two models do exist. For instance, in the discrete model, line solitons bifurcated from 
the X point in Fig. [5] are transversely stable for all ^ > iic ~ 21.5; while in the continuous model, soliton trains 
bifurcated from the X point in Fig. [2ja) can become unstable again when fj, gets close to the second Bloch band. The 
reason for this difference is that the discrete model can not capture the multi-band coupling, as we have mentioned 
earlier in this section. 



V. SUMMARY 



In summary, we have reported the existence of transversely stable soliton trains in optics. These soliton trains are 
found in two-dimensional square photonic lattices when they bifurcate from X points (with saddle-shaped diffraction) 
inside the first Bloch band, and their amplitudes are above a certain threshold. These stable soliton trains arise due 
to the combined effect of the photonic lattice, proper transverse phase relation, and strong longitudinal localization. 
We have also shown that soliton trains with low amplitudes or bifurcated from edges of the first Bloch band (F and 
M points) still suffer transverse instability. These results have been obtained in both the continuous lattice model 
and the discrete NLS model, and results from both models are in very good qualitative agreement. 
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